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Abstract. 

Due to the nonperturbative masslessness of the ghost field, ghost loops that contribute to 
gluon Green’s functions in the Landau gauge display infrared divergences, akin to those one 
would encounter in a conventional perturbative treatment. This is in sharp contrast with gluon 
loops, in which the perturbative divergences are tamed by the dynamical generation of a gluon 
mass acting as an effective infrared cutoff. In this paper, after reviewing the full nonperturbative 
origin of this divergence in the two-gluon sector, we discuss its implications for the three- 
and four-gluon sector, showing in particular that some of the form factors characterizing the 
corresponding Green’s functions are bound to diverge in the infrared. 


1. Introduction 

In the past few years, the infrared (IR) behavior of Yang-Mills Green’s functions in the Landau 
gauge has been the subject of numerous studies both in the continuum and on the lattice. 
Thanks to the use of a variety of different theoretical approaches (see, e.g., [T] and references 
therein) as well as extensive numerical simulations on large lattices O El 01 [B], a consistent 
picture has unequivocally emerged for the different n-point sectors of the theory. 

Most notably, it has been firmly established that the gluon propagator saturates at small 
momenta in a way consistent with the presence of a dynamically generated gluon mass PEIEIE]; 
the ghost propagator is instead essentially free in the same momentum region: in this case it is 
the ghost dressing function (defined as q'^ times the propagator, see below) that saturates to a 
finite non-vanishing value mm- 

This characteristic behavior of the two-point sector, which has been found to be valid in three 
and four space-time dimensions, for SU(3) and SU(2) gauge groups, and with or without the 
inclusion of dynamical quarks, profoundly affects the IR behavior of the theory’s entire tower of 
n-point Green’s functions [12]. In fact, it turns out that n-point functions exhibiting an ‘ancestor’ 
ghost-loop (i.e., a ghost-loop that is originally present in the lowest order perturbative expansion 
of the function under scrutiny), will develop a logarithmic IR singularity: the contribution of such 
diagram will correspond in fact to a pure logarithm, log (where ^ is some chosen scale), 

which is unprotected, in the sense that there is no mass term in its argument that could tame 
the corresponding divergence in the low momenta region. On the contrary, if the lowest order 
loops have a circulating gluon (and therefore a gluon propagator appears), the corresponding 


logarithm will be of the type log(g^ — thus, due to the presence of the dynamically 

generated gluon mass, such contributions will be finite for arbitrary low momenta. 

In the following, after reviewing the full nonperturbative analysis of the divergent ancestor 
ghost loop appearing in the 2-point sector and its effects on the gluon propagator and (inverse) 
dressing function, we will discuss the implications for the gluon three- and four-point sectors. In 
particular, we will show that the presence of these loops implies the appearance of IR divergences 
in some of the form factors comprising the general Lorentz decomposition of the corresponding 
vertices 

2. The 2-gluon sector 

The so-called PT-BFM framework [141 fT5l I16| . originating from the combination of the Pinch 
Technique (PT) [6l [T71 [181 [19] with the Background Field Method (BFM) [20], turns out to be 
particularly suited for studying the problem at hand, as it allows to separate in a gauge invariant 
way the ghost and gluon contributions to Green’s functions. Consequently, one is able to isolate 
ancestor ghost loops and study their IR behavior in a meaningful way. 

Within the PT-BFM framework, one considers the Schwinger-Dyson equation (SDE) of Fig.[Tj 
describing the propagator of a quantum (Q) and a background (B) gluon. Writing for 

the conventional (i.e., QQ) propagator 

~ )) ~ dtMu ~ Q.^iQi'/Q 1 ( 1 ) 

and similarly for A^,y(g), it turns out that the two propagators are related by a so-called 
Background-Quantum identity (BQI), reading (see again Fig. [T]) 

_ 6 

[1 + = A-^{q^)P^^{q) = q'^P^uiq) + i ( 2 ) 

i=l 

The auxiliary function G{(f‘) appearing above, corresponds to the metric form factor of a special 
Green’s function that is typical of this framework and describes the ghost-gluon dynamics. If 
one introduces the ghost dressing function F{q^) through 

F{q^) = q^D{q^), ( 3 ) 

D{q^) being the full ghost propagator, one has the approximate identity 122] 

l + G{q^)^F-\q^), ( 4 ) 

a relation that becomes exact at q^ = 0. 

The advantage of employing the BQI ([2]), and therefore considering the BQ self-energy 
diagrams rather than the QQ ones, resides in the fact that owing to the background Ward 
identity, all subsets of graphs enclosed within each box of Fig. [1] give rise to a transverse 
contribution [alia [IS]. Thus, their individual treatment, or, in fact, the total omission of 
entire subsets from one’s analysis, does not tamper with the transversality of the gluon self- 
energ}^. 

Now, in the presence, of an IR finite gluon propagator, the function A{q'^) can be decomposed 
as 

A-\q‘^)=q‘^J{q‘^)+m\q‘^), ( 5 ) 

^ In fact, it has been shown in |23| that the PT-BFM framework is not an ad-hoc procedure, rather it naturally 
emerges from the requirement of antiBRST invariance of the SU(N) Yang-Mills action. The additional identities 
(background Ward identity, BQIs, etc.) present in this novel formulation, are then none but manifestations of the 
underlying BRST-antiBRST invariance of the theory. 



Figure 1. (color online). The SDE obeyed by the QB gluon propagator. Each of the 
three different boxes (continuous, dashed, and dotted line) contains an individually transverse 
subgroup of diagrams, i.e., q^{{ai) + (aj+i)]^^. = 0 with i = 1, 2,3. Black (white) blobs represent 
fully dressed 1-PI (connected) Green’s functions; finally, small gray circles appearing on the 
external legs indicate background gluons. 


where is the inverse of the gluon dressing function, while is the dynamically 

generated gluon mass. Notice that Eq. ([2]) is satisfied independently by J{q^) and m?‘{q^), 
so that one has 

X(g2) = F{q^)X{q^)- X{q^) = J{q^),m\q^), (6) 

where we have used the approximation dH). 

The block-wise transversality property of the BQ propagator together with the BQI ([^, 
allows then for a meaningful separation of two kinds of contribution to J{q^)'. the ones stemming 
from the ghost graphs ( 03 ) and ( 04 ), and the ones stemming from the remaining gluon graphs. 
Denoting them respectively by Jc{q^) and Jg{q^), one has therefore 

q^Jg{q^)P^,u{q) = F{q^)[{ai) + {a 2 )]gu + F{q‘^)[{a 5 ) + {ae)]g^,, 

q^Jc{q^)Fg^{q) = E(g 2 )[(03)-p (04)1^1. ( 7 ) 

where on the right-hand side we assume that one is evaluating only the terms that vanishes as 
q^ goes to zero (the non vanishing terms contributing instead to the mass equation, see [ZlIB])- 
Then one finally has 

J(q^) = l + Jg{q^) + J,{q^), ( 8 ) 

with the “ 1 ” corresponding to the tree-level term. 

It turns out that there is a profound difference in the behavior of the gluon inverse dressing 
functions Jc{q^) and Jg{q‘^), which ultimately reflects the fact that the particle circulating in 
the loops of the corresponding diagrams have completely different behavior in the IR. In fact, 
in this momentum regime ghosts behave like massless free particles, D{q^) ~ F{^)/q^, whereas 
gluons are effectively massive, and therefore A(g^) ~ l/m^( 0 ). 

This difference between the inverse dressings Jc{q'^) and Jg{q'^) is most readily understood 
at the lowest order in the perturbative expansion: in the four-dimensional case, Jg{q^) develops 
a logarithm tamed by the presence of the effective IR cutoff provided by the dynamical gluon 
mass; on the contrary, Jc{q^) displays an unprotected logarithm, which vanishes at a finite value 
of then reverses its sign, becoming finally divergent at the origiij^. As a consequence, one can 

^ In the three-dimensional case the divergence in Jc{<l^) is linear in q with ~ 1/? while for the gluon 

Jg{q^) ~ arctan 5 / 2 m. 










easily show that the gluon propagator must display a maximum (and, consequently, its inverse 
a minimum) located in the deep IR region [12] . 

Consider now the full non perturbative case. The fact that the background gluon-ghost vertex 
appearing in the ghost block of Fig. [1] satisfies a QED-like Ward identity, furnishes a closed all¬ 
order expression for the longitudinal part of this vertex (which is not possible to obtain for the 
conventional gluon-ghost vertex). This leaves the transverse (automatically conserved) part of 
the vertex undetermined; however, under mild assumptions on the behavior of the form factor 
characterizing the latter vertex structure, one can show that the neglected term would give rise 
only to IR subleading contributions to Jc{<f) (see [T^j for details). 

Following the analysis presented in ca, one then finds that the ghost diagrams contribute 
to the gluon inverse dressing function the term 


where 
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( 10 ) 


In the equations above N represents the number of colors, and we have introduced the d- 
dimensional measure = //^/(27r)'^Jd'^/c, with q, the’t Hooft mass and e = 4 — d. Then, using 
the fact that T(0) = 0, in the deep IR region one finds the following behavior [T^j 
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+ 0{q^). 


( 11 ) 


The first integral appearing in the expansion of S{q^), contains the expected logarithmic 
divergence, as it can be easily seen by setting F{k‘^) = 1. Since the full F{k‘^) saturates at a 
constant value in the IR, its presence will not qualitatively modify the behavior of the integral; 
it will rather change its prefactor from 1 to F{0). The second integral in S{q‘^), and therefore 
T{q^), is subleading: introducing spherical coordinates and integrating by parts, makes manifest 
that this integral is simply proportional to F{0). In fact, Eq. ([8|) can be rewritten as 


J(g2) = J,V) + ^^V), 


( 12 ) 


where 


4iq^ 


g^N 

2{d-l) 



m 

k‘^{k + qY ’ 


(13) 


whereas J^^{q^) represents the IR subleading terms (including the terms generated by gluon 
graphs). 

The terms appearing in Eq. m can be next evaluated numerically by using as input a 
functional fit to the SU(2) |3] and SU(3) [I] quenched lattice data for the ghost dressing function. 
The results are shown in Fig. [21 where, as anticipated, the IR logarithmic divergence is clearly 
identified by the linear behavior (in log scale) of the S term above. 
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Figure 2. (color online). The ghost-loop contribution, (f‘Jc{<f‘)^ to the full gluon kinetic term 
q^J{q^) for the SU(2) (left) and SU(3) (right) gauge groups. 


To check whether or not the nonperturbative propagator has a maximum as a consequence of 
the divergence of its ghost contribution J^iq^) (US]), let us consider the derivative of the inverse 
gluon propagator that reads 

i^-Hq")]' = [q^Jiq")]' + m'{q^) = 4iq") + [r\q^) + q^J'iq'^)] - m'iq^). (14) 

Evidently, the quantity in brackets is subleading in the IR, while the fact that the propagator 
is decreasing in the UV, ensures that the above derivative is positive in this region. In addition, 
the dynamical equation governing m^(g^) is known [71[8l[9] and its solutions are monotonically 
decreasing and possess a finite derivative in the origin. Thus we conclude that the derivative 
above must reverse the sign at a point q^ where the propagator will display a maximum. 

Thus, according to the reasoning developed so far, and as a consequence of the masslessness 
of the ghost field, the lattice data for the gluon propagator ought to display a maximum, located 
in the (deep) IR region. As displayed in the insets appearing on the right panels of Fig. [3l such 
a maximum appears to be indeed encoded in the lattice data for A, which reveals a suppression 
of the deep IR points independently from the gauge group chosen. In those same insets we also 
plot different fitting curves in which the position of the maximum is varied. 

For reasons that will become clear in the next section, it is interesting to study the full kinetic 
term q^J{q^), an indirect knowledge of whicid can be acquired by evaluating the combination 

q^J{q^)=^-\q^)-m{q^), (15) 

where A(g^) is obtained from the aforementioned fits to the lattice while m‘^{q‘^) is obtained by 
solving the associated mass equation. Notice that one expects that the quantity (fT^ develops 
a minimum at a location qj, which in general however will not coincide with q^. 

The indirect determination of q^J{q^) from Eq. (lISp . using as basic input the family of 
curves for A(g^) obtained in the previous step is shown in Fig. [3j First we established that, 
when the latter curves are used as input to determine the solution to the mass equation, the 
resulting masses are completely independent of the location and the size of the maximum of the 

® The reason why we do not perform a direct SDE analysis of this quantity is because we do not have a satisfactory 
control over some of the basic ingredients appearing in the integral equation governing in particular, and 

most notably contrary to what happens for the Landau gauge mass equation, the equation for J(q^) involves the 
fully-dressed four-gluon vertex, whose structure is presently poorly known (see also Sect. I3.2ll . 
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Figure 3. (color online). The dynamical gluon mass (left panels) and the propagator’s full 
kinetic part A~^(g^) — m?‘{q^) (right panels) for the SU(2) (top) and SU(3) (bottom) gauge 
groups. Whereas the solutions of the mass equation are clearly insensitive to the presence of a 
maximum in the propagator, as shown here for two representative cases, the full kinetic term 
develops a negative minimum (^j), whose position is marked in the right panels by open (down) 
triangles. Insets show in all cases the IR behavior of the various propagator fits used as input, 
together with the corresponding lattice data of Refs. [21E] and m. 


propagator (left panels of Fig. [3]). Once the combination (fT^ is formed (right panels of Fig. [3]), 
we observe that the full kinetic term obtained vanishes at the origin, decreases in the deep IR, 
and reaches a negative minimum before crossing zero and turning positive (we mark for each 
curve the location of the corresponding minimum, qj). 

Summarising, the fact that the ghost field remains nonperturbatively massless, as opposed 
to the gluon which acquires a dynamically generated mass, implies unavoidably the existence 
of a negative IR divergence in the dimensionless co-factor J{q^) of the kinetic part of the gluon 
propagator. While this divergence, which originates exclusively from one-loop dressed diagrams 
involving a ghost loop, does not spoil the overall finiteness of the gluon two-point function, 
it affects it in two different ways: (i) it forces the appearance of a maximum in the gluon 
propagator located at q^^ and, correspondingly, (ii) it makes the full kinetic term develop a 
minimum at a location qj. 
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Figure 4. Ancestor ghost loops in the case of the three- and four-gluon sector. For later 
convenience we show the color and Lorentz indices of the external legs together with the 
momentum flow (all momenta are entering). 


3. The n-gluon sector 

The generality of the analysis performed in the previous section, suggests that the IR divergence 
appearing in the two-point sector, is likely to manifest itself in other Green’s functions that 
contain an ancestor ghost loop. Natural candidates are clearly the three- and four-point 
functions, due to the presence respectively of the triangular- and box-like diagrams shown 
in Fig. m 

In order to study this issue, in what follows we resort to the quantity usually employed at 
the non-perturbative level for studying n-point functions in the Landau gauge, namely the ratio 


(Hl ? * * * 1 Qn^ 


(51’ • • • > (qi) • • • {qi)Ki---Z (51, • • • , Qn) 

TZ'-'Z (51, • • • , (gi) • • • iqi)TZ:'Z (ft, • • • , 9n) ’ 


(16) 


What Eq. (fTUI) achieves is therefore the projection of the full vertex (gi, • • • ,g„) under 

scrutiny onto a particular tensor structure (o'!, • • • ,ft), factoring out at the same time 

external leg correction^. Such a ratio depends generally by the modulo of the n — 1 independent 
momenta and the (n — l)(n — 2)/2 angles between them. 


3.1. The three-gluon sector 

We start by considering the case of the three gluon vertex, defined according to (all momenta 
entering) 


^A-AbA-iPl,P2,P3) = -i9^fupiPl,P2,P3)-, 

^puf''iPl,P2,P3) = r^'^lg/iuiPl - P2)p + gup{P2 - P3)ii + 9pp{P3 - Pl)u], (17) 

where are the real and totally antisymmetric structure constants, satisfying the 

normalization condition _ ]\f^ab^ 

Choosing to project the vertex on its tree-level tensor structure, and considering for the 
momenta the so-called orthogonal configuration, corresponding to setting the angle between pi 
and p 2 to 7r/2 and then taking the limit p\ -A- 0, it can be shown that [12| 

Rf''\q\t),T^/2]=F{Q)[q^J{qZ]'+ R^ZgZ. ( 18 ) 

^ This definition allows, when data are available, for a direct comparison with lattice results, as in this case only 
connected (as opposed to 1-PI) Green’s functions can be measured. 
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Figure 5. (color online). Left panel: Prediction for the zero-crossing compared with the 
SU(2) ratio ,0,7r/2) measured on the lattice in four-dimensions [27]. The inset shows a 

logarithmic plot of the same quantity. Right panel: Same thing, but for the three-dimensional 
case EZj; notice that in this case the ghost divergence is linear in q, and therefore the lattice is 
able to resolve both the zero crossing and the negative divergence. 


where we have set = Pi, and the last term collects all subleading corrections not contained in 
the first one. Then the IR behavior of the ratio (11811 is driven once again by the gluon inverse 
dressing function so that 

Rf\q^0,7r/2] ~ F{0)4{q‘^). (19) 

Thus, the expectation is that R^^°^ in the orth^onal configuration will vanish at a point qo an 
then have a negative logarithmic IR divergencqj. Notice that an estimate for qo is provided by 
qj as the relation (11811 reveals; in particular, for the SU(2) gauge group the minimum of the full 
kinetic term provides the estimate qo ~ 45 MeV, while for SU(3) we obtain qo ~ 130 MeV. 

Our SU(2) result is compared with the behavior of Rp°\q‘^,0,7r/2) obtained from lattice 
simulations [27| on the left panel of Fig. [5j as one can see while there is indeed an indication 
that the zero crossing is going to happen, the actual value is located too deep in the IR to be 
resolved with current lattice volumes. On the other hand, we show on the right panel of the 
same figure the three-dimensional case, where the leading ghost divergence J44) is linear in 
momentum rather than logarithmic: in this case the zero crossing and divergent behavior are 
clearly resolved by the lattice and our prediction qo ~ 175 MeV compares reasonably well with 
the available data m- 

3.2. The four-gluon sector 

We next turn our attention to the four-gluon vertex which is defined according to (all momenta 
entering) 

^A-AiA-AiiPl,P2,P3,P4) = -ig'^'^^fupaiPU P2, PS, Pa)] 

+ f°'''''4^{gtia9up- gpuPpa)- ( 20 ) 


® For related studies on the three-gluon vertex, see |24l 1251126| . 













































Figure 6. The 18 diagrams contributing to the four-gluon vertex in the one-loop dressed 
approximation. The diagrams are divided in four different classes: (a) ghost boxes, (6) gluon 
boxes, (c) triangles, (d) fishnets (which carry a statistical factor of 1/2). 


As already remarked this is the most poorly understood vertex of the theory, e.g., no lattice 
simulation of this quantity exists to date (and, consequently, no data on any ratio in any 
momentum configuration are available). 

However, motivated by our successful description of the 2- and three-gluon sector, a 
preliminary nonperturbative study of this vertex can be attempted m- To this purpose, one can 
resort to a somewhat simplihed methodology, in which the different form factors are extracted 
directly from the evaluation of one-loop diagrams with fully dressed propagators but tree-level 
vertices (Fig. [6|). 

Even within this simplified setting, the calculation of the 18 one-loop dressed diagrams 
of Fig. [6] in a general momentum configuration would be a challenging task, due to the vast 
proliferation of tensorial structures. Indeed, at the level of rank-4 Minkowski tensors one has 
schematically the structures gg, gpiPj and PiPjPkPm, whereas for the rank-4 color tensors the 
possible combinations are of the type //, dd, fd and 65 {d being the totally symmetric SU(N) 
color tensor). This adds up to 138 possible tensors for a general kinematical configuration. 

Thus, in order to simplify the calculation as much as possible without compromising 
the physics we want to describe, we choose the momentum configuration {pi,P 2 ,P 3 ,P 4 :) = 
(P)P)P) “3p). This particular choice has the following advantages: 

(f) It gives rise to loop integrals that are symmetric under the crossing of external legs, thus 
reducing the amount of diagrams one needs to evaluate; 

(a) It allows to concentrate only on form factors multiplying tensor structures depending 
quadratically on the metric; all other possible structures will vanish when constructing 
the Landau gauge projectors R^; 

(Hi) It is the only momentum configuration in which the 1-PR contributions to the connected 
four-gluon Green’s function vanish, thus allowing to study the (projected) 1-PI component 
of the connected four-gluon vertex in isolatior]^. 

® This aspect would make the {p,p,p, —3p) as the configuration of choice in a possible attempt to study this 
vertex on the lattice. 


Within this configuration, for a general SU(N) gauge group one has 3 x 9 = 27 possible tensor 
structure^; however, for the special case of iV = 3, the additional identity [28] 


^ab^cd _|_ _|_ ^dd^bc _ ‘i\^clabrdcdr + dacrdbdr + dadrdbcr]} 


xacxbd 


xad xbc 


( 21 ) 


further reduces the total number of tensorial combinations down to 24. 

When all this is combined with the one-loop dressed approximation employed, it turns out 
that the tensor structures to be considered are in fact only two, as the result can be cast in the 
form 

= + vb(p 2 )G;“, (22) 

where is the tree-level terrsor defined in 1121 111 , while represents the totally symmetric 

tensor 

+ g,pgu. + gp.gup). (23) 

In addition, the IR leading term coming from the ghost diagrams (a) of Fig. ?? only contributes 
to the latter structure, as one has 


^abcd 

fii/pa 


{p,p,p,-3p) 
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i=l 
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A{p^) = 

Notice that as p^ —7 
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2 d 2 _ 1 


[ k^ 

■ {k-p)'^' 
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J^'ApZ 


0, Eq. ([2^ yields 
A{p^^ 


F{k)F{k +p)F{k + 2p)F{k + 3p) 
(k + p)"^ {k + 2py {k + 3p)^ 

[F^{k‘^) 


p2^o 2d{d + 2) J). k^ 


(24) 


(25) 


(26) 


that is, the form factor Vg{p^) will develop a logarithmic IR divergence (in four-dimensions). 
The (rather lengthy) expressions for the remaining class of diagrams can be found in [13j . 

The different contributions to the four-gluon vertex form factors can be evaluated using the 
functional fits to the quenched lattice data for the gluon and ghost two-point function^. The 
results are shown in Fig. [71 

Consider first the left panel, where we show the contributions to Ip(o). These are purely 
gluonic in nature, and all of them saturate in the IR. In particular we notice that, the contribution 
of the gluon boxes is negligible (indeed, as p ^ 0 it vanishes); the triangle terms feature a bump 
in the low momentum region, while the fishnet is negative. Adding everything up, one obtains 
the shape shown by the black line. 

In the case of Vq the situation is completely different (Fig. [71 right panel). Gluon contributions 
are again saturating in the IR; however, in this case, the ghost boxes take over below few 
hundreds MeV^, driving Vq to an IR logarithmic divergence. As far as the remaining diagrams 
are concerned, gluon boxes are negative in this case; in addition, they are almost perfectly 
cancelled by the triangle contributions. When the negative contribution from the fishnet 
diagrams is finally added, one obtains the shape shown by the black line of Fig. [71 in particular, 

^ The terms quadratic in the metric gives the 3 possible combinations g/ivgptr, g^pgua, and gpagvp', for the color 
structures on has 15 possibilities and 6 identities [28], and therefore 9 independent tensors. 

® In the case of the gluon propagator we have also employed a fit featuring the IR maximum discussed in Sect. [2] 
the results obtained were, however, independent from its presence. 













Figure 7. (color online). Individual one-loop dressed contributions to the tensor structure 
(left panel) and (right panel). The black line coincides with the coefficient Vp(o) (p^) 

and Vg{p^) respectively. 


notice the presence of a zero crossing, a feature that was also present in the ^^3 ratio previously 
studied in the case of the three-gluon vertex. 

Now, the result ([2^ is due to the one-loop dressed approximation employed; a general analysis 
based on Bose symmetry m, shows that the terms quadratic in the metric contributing to the 
four-gluon vertex in the {p,p,p, —3p) momentum configuration allow for an extra tensor, namely 
one has 


^abcd 

fiupa 


{p,p,p, -3p) 


■r r / 


(27) 


wher^ 




(28) 


This means that a possible lattice evaluation of the connected four-gluon function in this 
particular momentum configuration completely determines the terms of the four-gluon vertex 
quadratic in the metric tensor, through the measurement of the ratios 

Rf\p^) = F^(o)(p2) + lv3,,(/); Rf(p^) = Fc(/); (p2) = Vx'ip^) + -^V^,o) (p^). 

(29) 

According to our description one expects 17 ( 0 ) (p^) to be finite and Vg{p^) to display an IR 
divergence; nothing can be however said at the moment for the form factor Vxb however 
the vanishing of this latter quantity in the one-loop dressed approximation, points towards 
its finiteness. Thus, we would expect the measurement of only one divergent ratio, and namely 

Rf(p 2 ). 

® Bose symmetry alone does not permit to fix completely the tensor its explicit form 1281) is obtained by 

requiring that the latter should be orthogonal to the tensor G, that is 
























4. Conclusions 

The gluon and ghost field display a very different behavior in the deep IR: the latter remains 
nonperturbatively massless, whereas the former acquires a dynamically generated mass. This 
fact, which has been unequivocally established in the Landau gauge using discrete as well as 
continuous methods, turns out to have a profound impact on the n-gluon sector of the theory, 
as diagrams involving ghost loops gives unavoidably origin to IR divergences. 

In the case of the 2-gluon sector it is the gluon inverse dressing function J{q^) that shows such 
a divergence (with J{q^) ~ logg^ in four dimensions); while the presence of such a divergence 
does not interfere with the finiteness of the gluon 2 -point function (for the gluon full kinetic 
term is multiplied by q^), it nevertheless implies that the full propagator has an IR maximum 
located at g = and, correspondingly, q^J{q^) has a minimum, located at g = qj. 

For the three-gluon sector a (negative) IR divergence appears when projecting the full 
three-gluon vertex onto its tree level value, and choose the so-called orthogonal momentum 
configuration. Due to the relation (fT^ . the location of the point g = go at which ^ 3 ^°^ (^^,0, 7 r/ 2 ) 
crosses zero and turns negative can be roughly estimated from qj and turns out to be around 
130 MeV (45 MeV) for the N = 3 {N = 2) case. 

An IR divergence appears also when evaluating the four-gluon vertex in the {p,p,p, —3p) 
momentum configuration, even though in this case it does not manifest in the projector (p^) 
onto the tree-level tensor, rather in Rf{p^), where G is the totally symmetric tensor (I23p . 

The picture presented here, elaborated within the PT-BFM formalism which allows for a 
gauge-invariant separation of ghost and gluon contributions to the gluon propagator, is found 
to be in agreement with lattice data whenever the latter are available. 

Two are the questions that needs to be addressed. 

To begin with, since Green’s functions depend on the gauge fixing employed, it would be 
important to evaluate them in different gauges in order to ascertain what aspects of their 
nonperturbative behavior are affected by a change of gauge. This is particularly relevant in the 
2-point sector, as a recent preliminary study using a combination of SDE and Nielsen identitieJ^ 
has revealed that in the renormalizable ^ gauges the ghost dressing function vanishes in the deep 
IR [3l]. If this result persists refined studies, possibly including lattice simulations, its impact 
on the IR behavior of ancestor ghost loops (and consequently the n-gluon sector of the theory) 
needs to be thoroughly addressed. 

Second, one would like to see wether or not the presence of zero crossings and IR divergences 
has some impact on hadron phenomenology {e.g., the hadron spectrum), in order to connect 
results on the theory’s most basic building blocks with its observables properties, along the lines 
recently discussed in [32]. At a first sight it would look like the zero crossing momentum go is 
located too deep in the IR for both the three- and four-gluon vertex to affect the Bethe-Salpeter 
equations one needs to solve. However, a preliminary analysis of unquenching effects shows 
that dynamical quarks have the tendency to move go closer to the phenomenologically relevant 
region of few hundreds MeV. This is relatively easy to understand, as the main effect of adding 
dynamical quarks is to suppress the saturation point of the gluon propagator while leaving the 
ghost dressing function practically invariant laiMlIMI; consequently gluonic contributions will 
be suppressed whereas ghost contributions will be of the same size of the unquenched ones, 
which results in pushing the zero crossing towards higher momentum values with respect to the 
quenched case. 

These identities express the gauge-dependence of ordinary Greens functions (propagators, vertices, etc.) in 
terms of special auxiliary functions associated with an extended BRST sector |29[ l,S0 |. 
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